ABSTRACT
A Compton camera, which is based on the geometrical interpretation of Compton scattering, is a very promising gamma-ray imaging device considering its several advantages over the conventional gamma-ray imaging devices: high imaging sensitivity, 3-D imaging capability from a fixed position, multi-tracing functionality, and almost no limitation in photon energy. In the present study, a Monte Carlo-based, user-friendly Compton imaging simulator was developed in the form of a graphical user interface (GUI) based on Geant4 and MATLAB (TM). The simulator was tested against the experimental result of the double-scattering Compton camera, which is under development at Hanyang University in Korea. The imaging resolution of the simulated Compton image well agreed with that of the measured image. The imaging sensitivity of the measured data was 2~3 times higher than that of the simulated data, which is due to the fact that the measured data contains the random coincidence events. The performance of a stacking-structure type Compton camera was evaluated by using the simulator. The result shows that the Compton camera shows its highest performance when it uses 4 layers of scatterer detectors.
Subject(s)
Imaging, Three-Dimensional , KoreaABSTRACT
PURPOSE: In this study we propose a block-iterative method for reconstructing Compton scattered data. This study shows that the well-known expectation maximization (EM) approach along with its accelerated version based on the ordered subsets principle can be applied to the problem of image reconstruction for Compton camera. This study also compares several methods of constructing subsets for optimal performance of our algorithms. MATERIALS AND METHODS: Three reconstruction algorithms were implemented; simple backprojection (SBP), EM, and ordered subset EM (OSEM). For OSEM, the projection data were grouped into subsets in a predefined order. Three different schemes for choosing nonoverlapping subsets were considered; scatter angle-based subsets, detector position-based subsets, and both scatter angle- and detector position-based subsets. EM and OSEM with 16 subsets were performed with 64 and 4 iterations, respectively. The performance of each algorithm was evaluated in terms of computation time and normalized mean-squared error. RESULTS: Both EM and OSEM clearly outperformed SBP in all aspects of accuracy. The OSEM with 16 subsets and 4 iterations, which is equivalent to the standard EM with 64 iterations, was approximately 14 times faster in computation time than the standard EM. In OSEM, all of the three schemes for choosing subsets yielded similar results in computation time as well as normalized mean-squared error. CONCLUSION: Our results show that the OSEM algorithm, which have proven useful in emission tomography, can also be applied to the problem of image reconstruction for Compton camera. With properly chosen subset construction methods and moderate numbers of subsets, our OSEM algorithm significantly improves the computational efficiency while keeping the original quality of the standard EM reconstruction. The OSEM algorithm with scatter angle- and detector position-based subsets is most available.